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We suggest that non-trivial correlations between the dark matter particle mass and collider based 
probes of missing transverse energy H^^^^ rnay facilitate a two tiered approach to the initial discovery 
of supersymmetry and the subsequent reconstruction of the lightest supersymmetric particle (LSP) 
mass at the LHC. These correlations are demonstrated via extensive Monte Carlo simulation of 
seventeen benchmark models, each sampled at five distinct LHC center-of-mass beam energies, 
spanning the parameter space of No-Scale J-"-SC/(5). This construction is defined in turn by the 
union of the J^-lipped 5'(7(5) Grand Unified Theory, two pairs of hypothetical TeV scale vector- 
like supersymmetric multiplets with origins in J^-theory, and the dynamically established boundary 
conditions of No-Scale Supergravity. In addition, we consider a control sample comprised of a 
standard minimal Supergravity benchmark point. Led by a striking similarity between the i/™'"" 
distribution and the familiar power spectrum of a black body radiator at various temperatures, 
we implement a broad empirical fit of our simulation against a Poisson distribution ansatz. We 
advance the resulting fit as a theoretical blueprint for deducing the mass of the LSP, utilizing only 
the missing transverse energy in a statistical sampling of > 9 jet events. Cumulative uncertainties 
central to the method subsist at a satisfactory 12-15% level. The fact that supersymmetric particle 
spectrum of No-Scale J--SU{5) has thrived the withering onslaught of early LHC data that is steadily 
decimating the Constrained Minimal Supersymmetric Standard Model and minimal Supergravity 
parameter spaces is a prime motivation for augmenting more conventional LSP search methodologies 
with the presently proposed alternative. 

PACS numbers: ll.lO.Kk, ll.25.Mj, ll.25.-w, 12.60.Jv 



I. INTRODUCTION 



A. Dark Matter and Missing Transverse Energy 



High-precision experiments have revealed that the 
matter content of our Universe is predominantly non- 
baryonic and dark. The presence of large scale structure 
formation indicates that dark matter should be primar- 
ily cold, i.e. possessing a non-relativistic thermal speed 
at the time of its own decoupling. Measurements of 
the anisotropics in the Cosmic Microwave Background 
(CMB) radiation, most recently and conclusively by the 
WMAP [1H3 satellite experiment, allocate the propor- 
tional energy densities of baryonic matter, cold dark mat- 
ter (CDM), and dark energy (a cosmological constant) at 
about 4%, 23%, and 73% respectively. A stable Weakly 
Interacting Massive Particle (WIMP) with a canonical 
weak-scale mass can naturally have a thermal CDM relic 
density at the correct magnitude order [4] . Although the 
Standard Model (SM) does not possess a CDM candi- 
date, suitable particles are manifest in various SM ex- 
tensions, most notably in constructions featuring super- 
symmetry (SUSY), which itself also efhciently addresses 
the quantum stability of the electroweak (EW) gauge hi- 



erarchy. Decay of the Lightest Supersymmetric Particle 
(LSP), typically associated with a neutralino eigenstate, 
may be blocked by a discrete i?-parity symmetry, result- 
ing in an exceedingly natural CDM candidate [1, Q . Al- 
ternative examples such as the little Higgs model with 
T-parity [y, 0] and extra dimensional models with KK 
parity is] are not considered further here. 

In general, if dark matter particles are produced in a 
collider, being [/(1)em neutral, SU{3)c singlets, they will 
escape the detector without direct observation. The sole 
evidence of their production and rapid departure will be 
that of a deficit in the zero sum of momentum over all 
observed particles with a component trajectory which is 
perpendicular to the beamline. This so-called "missing 
transverse energy" (H^^^^) is measurable to fine precision 
at the LHC, which has been steadily accumulating data, 
thus far without any positively reported evidence for su- 
persymmetry, from y/s = 7 TeV proton-proton collisions 
since March 2010 [MtIII- In the current early opera- 
tional stage of the LHC, our attention has unquestionably 
been piqued by the relentless severe reductions inflicted 
upon the viable minimal Supergravity (mSUGRA) and 
Constrained Minimal Supersymmetric Standard Model 
(CMSSM) parameter space. It was concluded by one 
particular study that the miniscule 35 pb~ of initially 



reported luminosity was alone sufBcient to impeach more 
than 99% of the parameter space [12|1 , although other an- 
alysts remained more circumspect [1^. The negative re- 
sults have now held through the collection of more than 
a femtobarn of integrated luminosity by each of the AT- 
LAS 14] and CMS [lH collaborations, although we must 
remark that tantalizing low-statistics excesses do exist 
in both cases, which are neatly accounted for by simula- 
tion of the No-Scale T-SU{5) collider-detector signal [l6| . 
Regardless, should this situation much longer persist into 
the deconstruction of the burgeoning wealth of collected 
raw data, reported to have already reached the impres- 
sive milestone of 5 fb~ by the close of 2011, then the first 
great question of the LHC era shall become significantly 
heightened. Namely, that is the question of whether there 
exist specifically detailed SUSY and/or superstring based 
post-SM constructions which simultaneously escape the 
existing experimental onslaught, yet also remain in prin- 
ciple discoverable to the present generation collider. 

Granting, momentarily, an answer in the tacit affirma- 
tive to the prior question of discovery, one may speculate 
regarding the equally great question immediately subse- 
quent, that being the issue of the detailed SUSY CDM 
classification. Differentiation between contending scenar- 
ios, both within the CMSSM and without, will require a 
probe at present and future collider and direct detection 
experimentSjthe latter now led by the XENONIOO col- 
laboration [13], of the intrinsic properties of the CDM, 
such as its spin and mass, and the structure of its in- 
teraction with SM particles. Assuming a SUSY LSP 
to represent the CDM, we may in principle make a col- 
lider based measurement of the LSP mass from the decay 
chains of heavier particles, for example linking the squark 
and gluino to the lightest neutralino in the supersym- 
metric SM [l^, [I^l , though this method is highly model 
dependent, and there are myriad complications with its 
interpretation. Firstly, since the momentum component 
along beamline can by its nature never be collected, it 
is impossible to reconstruct the true net missing energy 
total. Secondly, i?-parity implies that the LSP candi- 
date will always ultimately appear pairwise. Moreover, 
the two neutralinos may be kinematically decoupled from 
the original back-to-back pair production event, e.g. of 
two heavy squarks qq, by an intricate decay cascade, 
such that the final relative angular orientation is purely 
random, allowing for any possible partial cancellation. 
Finally, the dark matter component of the total miss- 
ing energy H^^^^ must be segregated from other poten- 
tial sources, such as energetic neutrinos and energy mis- 
measurement, for example by particle escape through a 
calorimeter seam. Although the prospect of a substan- 
tial loss of the latter type is unlikely, it competes with 
the ever more extraordinarily rare target events. These 
foreboding circumstances compel some candid apprehen- 
sion against the hope that any fundamental correlation 
might endure between the dark matter particle mass and 
the missing transverse energy measurement. 

In the present work we set out to directly confront, 



by the aid of a massive Monte Carlo collider and de- 
tector simulation, the two great questions here posed. 
The two tiered nature of our presentation, approach- 
ing first the prospect of SUSY CDM discovery, and sec- 
ondly the subsequent classification of the SUSY CDM 
signal via explicit numerical correlation of the LSP mass 
with the peak of the missing transverse energy distri- 
bution, is in correspondence with our view of how the 
experimental reality might plausibly unfold in the com- 
ing months and years. Our study takes the context of 
a specific model, dubbed No-Scale J'-SUib) ^M^, 
representing the merger of the J^- lipped SU{5) Grand 
Unified Theory (GUT) (H-Ill], two pairs of hypotheti- 
cal TeV scale vector-like supersymmetric multiplets with 
origins in J^-theory |35l439l |. and the dynamically estab- 
lished boundary conditions of No-Scale Supergravity [40!- 
|44| . The construction features an essential phenomeno- 
logical consistency, profound predictive capacity, a sin- 
gularly distinctive experimental signature, and imminent 
testability. The core of the model space remains exper- 
imentally unblemished, in stark contrast to the incisive 
cuts already made into the heart of the more traditional 
mSUGRA/CMSSM formulations. This deft evasion of 
all existing LEP, Tevatron, and LHC sparticle mass lim- 
its stems from the characteristic sparticle mass hierarchy 
TOj < rrig < rriq of a light stop and gluino, both com- 
fortably lighter than all other squarks, a feature which 
is stable across the full model space. If the tide advanc- 
ing against the CMSSM should continue unabated, then 
island refuges such as No-Scale J^-SU{5) will demand in- 
creasing deference in the search for a physically feasible 
supersymmetric grand unified theory. 



B. A No-Scale T-SUi5) Primer 

The No-Scale J^-SU{5) construction inherits all of the 
most beneficial phenomenology [45] of flipped SU{5) |32| - 
|34| . including fundamental GUT scale Higgs represen- 
tations (not adjoints), natural doublet-triplet splitting, 
suppression of dimension-five proton decay and a two- 
step see-saw mechanism for neutrino masses, as well as 
all of the most beneficial theoretical motivation of No- 
Scale Supergravity |40l - l44| . including a deep connection 
to the string theory infrared limit (via compactification of 
the weakly coupled heterotic theory [4g or M-theory on 
S^ jZi at the leading order [431), the natural incorpora- 
tion of general coordinate invariance (general relativity), 
a mechanism for SUSY breaking which preserves a van- 
ishing cosmological constant at the tree level (facilitating 
the observed longevity and cosmological flatness of our 
Universe [40|), natural suppression of CP violation and 
flavor-changing neutral currents, dynamic stabilization of 
the compactifled spacetime by minimization of the loop- 
corrected scalar potential and a dramatic reduction in 
parameterization freedom. 

Written in full, the gauge group of flipped S'C/(5) is 
SU(h) X U{\)x, which can be embedded into S'O(IO). 



The generator U{1)y' is defined for fundamental five- 
plets as —1/3 for the triplet members, and +1/2 for 
the doublet. The hypercharge is given by Qy = (Qx — 
Qy')/5. There are three families of Standard Model (SM) 
fermions, whose quantum numbers under the SU{5) x 
U{l)x gauge group are 

f, =.(10,1) ; /, -(5,-3) ; T, = (1,5), (1) 

where i — 1,2,3. There is a pair of ten-plet Higgs for 
breaking the GUT symmetry, and a pair of five-plet Higgs 
for electroweak symmetry breaking (EWSB). 



if = (10,1) 

h^ (5,-2) 



H = (10,-1) 

^=(5,2) 



(2) 



Since we do not observe mass degenerate superpartners 
for the known SM fields, SUSY must itself be broken 
around the TeV scale. In the minimal supergravities 
(mSUGRA), this occurs first in a hidden sector, and 
the secondary propagation by gravitational interactions 
into the observable sector is parameterized by universal 
SUSY-breaking "soft terms" which include the gaugino 
mass Ml/2, scalar mass Mq and the trilinear coupling 
A. The ratio of the low energy Higgs vacuum expec- 
tation values (VEVs) tan/3, and the sign of the SUSY- 
preserving Higgs bilinear mass term /i are also undeter- 
mined, while the magnitude of the fj, term and its bilin- 
ear soft term _B^ are determined by the Z-boson mass 
Mz and tan/3 after EWSB. In the simplest No-Scale sce- 
nario, Mo=A=_B^=0 at the unification boundary, while 
the complete collection of low energy SUSY breaking 
soft-terms evolve down with a single non-zero parameter 
Ml/2. Consequently, the particle spectrum will be pro- 
portional to Ml/2 at leading order, rendering the bulk 
"internal" physical properties invariant under an overall 
rescaling. The rescaling symmetry can likewise be bro- 
ken to a certain degree by the vector-like mass parameter 
My, although this effect is weak. 

The matching condition between the low-energy value 
of Bf^ that is demanded by EWSB and the high-energy 
B^ —Q boundary is notoriously difficult to reconcile un- 
der the renormalization group equation (RGE) running. 
The present solution relies on modifications to the /3- 
function coefficients that are generated by the inclusion 
of the extra vector-like multiplets, which may actively 
participate in radiative loops above their characteristic 
mass threshold Afy- Naturalness in view of the gauge 
hierarchy and /i problems suggests that the mass My 
should be of the TeV order. Avoiding a Landau pole for 
the strong coupling constant restricts the set of vector- 
like multiplets which may be given a mass in this range to 
only two constructions with flipped charge assignments, 
which have been explicitly realized in the i^-theory model 
building context [35l - l37| . We adopt the multiplets 

XF = (10, 1) = (Xg, XD", XN") ; XF = (TO, -1) 

XI = (1,-5) ; XI=(1,5)=XS^ (3) 



where XQ, XD", XE'^ and XN'^ carry the same quan- 
tum numbers as the quark doublet, right-handed down- 
type quark, charged lepton and neutrino, respectively. 
Alternatively, the pair of SU{5) singlets may be dis- 
carded, but phenomenological consistency then requires 
the substantial application of unspecified GUT thresh- 
olds. In either case, the (formerly negative) one-loop 
/3-function coefficient of the strong coupling as becomes 
precisely zero, flattening the RGE running, and generat- 
ing a wide gap between the large 032 — 0:3 (Mz) ~ 0.11 
and the much smaller ax at the scale M32 of the inter- 
mediate flipped SU{5) unification of the SU{3) x 5'C/(2)l 
subgroup. This facilitates a very significant secondary 
running phase up to the final SU{5) x U{l)x unifi- 
cation scale [38], which may be elevated by 2-3 or- 
ders of magnitude into adjacency with the Planck mass, 
where the B., = boundary condition fits like hand to 
glove [13, EE [H. This natural resolution of the "little 
hierarchy" problem corresponds also to true string-scale 
gauge coupling unification in the free fermionic string 
models 1351 . ISOl or the decoupling scenario in F-theory 
models [3y, [33| , and also helps to address the monopole 
problem via hybrid inflation. 

The modifications to the /3-function coefficients from 
introduction of the vector-like multiplets have a paral- 
lel effect on the RGEs of the gauginos. In particular, 
the color-charged gaugino mass M3 likewise runs down 
fiat from the high energy boundary, obeying the rela- 
tion M3/M1/2 ~ a3(Mz)/a3(M32) ~ 0(1), which pre- 
cipitates a conspicuously light gluino mass assignment. 
The S'C/(2)l and hypercharge U{1)y associated gaug- 
ino masses are by contrast driven downward from the 
Ml/2 boundary value by roughly the ratio of their cor- 
responding gauge couplings (a2,aY) to the strong cou- 
pling as- The large mass splitting expected from the 
heaviness of the top quark via its strong coupling to the 
Higgs (which is also key to generating an appreciable ra- 
diative Higgs mass shift A irif^ [3l|) is responsible for 
a rather light stop squark ti. The distinctively predic- 
tive TOjj^ < rrig < rriq mass hierarchy of a light stop and 
gluino, both much lighter than all other squarks, is stable 
across the full No-Scale J^-SU{5) model space, but is not 
precisely replicated in any phenomenologically favored 
constrained MSSM (CMSSM) constructions of which we 
are aware. 

This spectrum generates a unique event topology start- 
ing from the pair production of heavy squarks qq, ex- 
cept for the light stop, in the initial hard scattering pro- 
cess, with each squark likely to yield a quark-gluino pair 
q — >■ qg. Each gluino may be expected to produce events 
with a high multiplicity of virtual stops, via the (possibly 
off-shell) g — 7> t transition, which in turn may terminate 
into hard scattering products such as -^ W~^W'bbxi 
and W~bbT'^i'rXiJ where the W bosons will produce 
mostly hadronic jets and some leptons. The model de- 
scribed may then consistently exhibit a net product of 
eight or more hard jets emergent from a single squark 
pair production event, passing through a single inter- 



mediate gluino pair, resulting after fragmentation in a 
spectacular signal of ultra-high multiplicity final state 
jet events. We remark also that the entirety of the vi- 
able T-SUlSi) parameter space naturally features a dom- 
inantly Bino LSP, at a purity greater than 99.7%, as is 
exceedingly suitable for direct detection, for example by 
XENONIOO [13, HI- There exists no direct Bino to Wino 
mass mixing term. This distinctive and desirable model 
characteristic is guaranteed by the relative heaviness of 
the Higgs bilinear mass /x, which in the present construc- 
tion generically traces the universal gaugino mass M1/2 
at the boundary scale Mjr, and subsequently transmutes 
under the RGEs to a somewhat larger value at the elec- 
troweak scale. 

A majority of the bare-minimally constrained 26] pa- 
rameter space of No-Scale J-'-SU{5), as defined by con- 
sistency with the world average top-quark mass rrit, the 
No-Scale boundary conditions, radiative EWSB, the cen- 
trally observed WMAP7 CDM relic density j3], and pre- 
cision LEP constraints on the lightest CP-even Higgs bo- 
son rrih [5l|, [52| and other light SUSY chargino and neu- 
tralino mass content, remains viable even after careful 
comparison against the first inverse femtobarn of LHC 
data (Tg. l29ll . Moreover, a highly favorable "golden" sub- 
space [20, UM, [SJl exists which may simultaneously ac- 
count for the key rare process limits on the muon anoma- 
lous magnetic moment {g — 2)^ and the branching ra- 
tio of the flavor-changing neutral current decays b ^ s"f 
and -Bg — > fj,'^fj,~. The intersection of these experimen- 
tal bounds is highly non-trivial, as the tight theoreti- 
cal constraints, most notably the vanishing of B^ at the 
high scale boundary, render the residual parameteriza- 
tion deeply insufficient for arbitrary tuning of even iso- 
lated predictions, let alone the union of all predictions. In 
addition, a top-down consistency condition on the gaug- 
ino boundary mass Afx/2, and the parametrically cou- 
pled value of tan/3, is dynamically determined at a sec- 
ondary local minimization dVYnin/dMi^2 = of the mini- 
mum of the Higgs potential Kiiin; since M1/2 is related to 
the modulus field of the internal string theoretic space, 
this also represents a dynamic stabilization of the com- 
pactification modulus. The result is demonstrably con- 
sistent with the bottom-up phenomenological approach 
[23 . I23I . |26| . |30 | , and we note in particular a rather dis- 
tinctive conclusion that is enforced in both perspectives: 
the ratio tan/3 must have a value very close to 20. 



II. TIER 1: SUSY CDM DISCOVERY 

A. Missing Transverse Energy as a SUSY Signal 

In the first tier of this study, we present an exhaustive 
Monte Carlo simulation of seventeen benchmark points, 
as detailed in Table HI Three heavier benchmarks have 
been chosen here to supplement the original fourteen 
from Ref. \26\. Together, these sample points span the 
parameter space of bare minimal phenomenological con- 



straints [201 on No-Scale T-SU{5). We also consider the 
third "Snowmass Points and Slopes" (SPS) benchmark 
point (SPS3) [Hi as a well established mSUGRA con- 
trol sample, with tolsp = 161.7 GeV, Mo = 90 GeV, 
A/1/2 == 400 GeV, TOt = 175 GeV and tan/3 = 10. 
We present an example SUSY spectrum for the favor- 
able mLSP = 100.0 GeV, M1/2 = 520 GeV benchmark in 
Table [TTj 

A major new feature of the current study is the exten- 
sion of our prior simulation coverage to include five dif- 
ferent LHG center-of-mass beam energies, starting from 
the presently operational y/s = 7 TeV machine, and cul- 
minating in the bright future of a ^/s = 14 TeV LHC. 
Our initial purpose here is the development of a compre- 
hensive map of the prospects for SUSY discovery across 
this model and energy space, based on the raw count 
of events surviving a set of cuts designed for suppres- 
sion of SM backgrounds. The entirety of our collider 
based analysis is in fact dependent upon adoption of an 
ultra-hi gh m ultiplicity (> 9) jet event cutting methodol- 
ogy /S, [23 • The heart of this approach is the observa- 
tion that the decay chains described in the introduction, 
which are prevalent for the signature T-SU{5) mass hi- 
erarchy, will very frequently result in at least eight hard 
jets per SUSY event (the W bosons will produce mostly 
hadronic jets and some leptons), and given fragmenta- 
tion processes, perhaps a much larger number of softer 
final state jets. Since the ultra-high jet regime has been 
demonstrated [2a| to be naturally suppressed in the SM 
backgrounds, we are further able to relax certain of the 
harsh cuts which are very effective for separating out the 
CMSSM in intermediate jet searches, but which simulta- 
neously exert a costly attrition against our signal. 

There are two varieties of related selection cuts which 
we will entertain. The original ULTRA selections, so 
named for their emphasis of the ultra-high jet multiplic- 
ity signal, are carefully documented in Ref. [25 • These 
cuts are characterized by a rather low transverse mo- 
mentum threshold of px ^ 20 GeV per jet, and require 
at least 150 GeV of missing energy H^^^^. Any cut on 
the CMS collaboration statistic ax, which is designed 
to discriminate against false missing energy signals, is 
expressly forgone, as we have found it to be likewise sys- 
tematically biased against the ultra-high jet count signal. 
Again, the high jet count threshold is in itself a quite 
effective discriminant against the SM background. How- 
ever, the ULTRA selection cuts have a rather substantial 
drawback; specifically, we have no deeply reliable esti- 
mates of the precise amount of background signal which 
does penetrate. Although we have a baseline established 
by substantial Monte Carlo of the crucial ti + jets SM 
component, the massive cross-section of the QCD con- 
tribution makes a complete proprietary simulation pro- 
hibitively difficult. 

We thus also consider a second variety of selection 
cuts, which we refer to as the CMS HYBRID, that has the 
distinct advantage of pre-tabulated SM background es- 
timates from the CMS collaboration itself. We describe 



TABLE I: We tabulate various statistics for seventeen representative points selected from the viable No-Scale J^-SU{5) param- 
eter space, each satisfying the bare-minimal phenomenological constraints outlined in Ref. 26]. The dimonsionful parameters 
JT^LSP, A4"i/2, Afv,nit, _B^(Afjr), ^(Afz), and mo are in units of GeV. The discoverability statistic S/y/B + 1 is computed at a 

y^ — 7 TeV LHC, scaled to 1 fb^^ of integrated luminosity, employing the CMS HYBRID selection cut methodology 16] in our 
own Monte Carlo, and adopting a corresponding SM background of B = 3.1 events, as extracted from a report by the CMS 
collaboration 1 1511. 
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445 


1125 


19.9 


174.4 


0.58 


789 
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3700 
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20.7 
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-0.19 
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1.60 
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0.94 
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0.1131 
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108.8 


1.4 
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0.96 


133.2 
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22.5 
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0.46 


1143 


0.1118 


0.59 


155.9 


0.1 
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6000 


23.0 
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0.63 


1335 


0.1120 


0.33 


190.5 


0.0 



TABLE II: The detailed SUSY particle spectrum (in GeV) 
for the experimentally favorable benchmark with M1/2 = 520 
GeV, tan^g = 20.7, Mv = 1725, and mt = 174.4 GeV. Note 
that the presented light CP-even Higgs mass rrih does not 
include radiative corrections from the vector-like "flippon" 
multiplets, which may account for an upward shift of around 
3-4 GeV ^. 
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the details of this selection methodology in Refs. |16l.l29l|. 
where it is compared to the CMS report [l^ • By carefully 
extracting just those events above the nine jet threshold 
from the detailed histogram plot presented in that ref- 
erence, we estimate a SM background of approximately 
3.1/fb~^ events for the CMS HYBRID selection cuts. In 
most critical regards, and essentially in suppression of 
the ax selection, these cuts are quite similar to our orig- 
inal ULTRA preference, which is the reason for the "hy- 
brid" description. There are various incidental historical 
distinctions, including a reduction of the missing energy 
threshold to H:^'''^ > 100 GeV, and the introduction of 
cuts on the electromagnetic fraction and an upper bound 
on the ratio of missing energy carried by hard to soft jets. 
The most notable difference is an increase in the trans- 



verse momentum required per hard jet, to pT > 50 GeV. 
While the greater overall jet count produced by the ULTRA 
cuts is quite valuable for discovery at the lighter end of 
the J^-SU{5) model space (well below M1/2 ~ 500 GeV 
or tolsp ~ 100 GeV), we have found [ig] that the greater 
background suppression of the CMS HYBRID cuts becomes 
valuable for the lower-producing heavier spectra. Since 
our presently favored [T^ region of the model is near the 
just mentioned transition, where the discoverability pro- 
duced by the two cuts appears to come roughly into par- 
ity, we consider the CMS HYBRID selections to be a very 
acceptable alternative. Because of the clear advantage in 
quantification of the background competition, we opt to 
employ these cuts in the first tier of our study, dealing 
with model discoverability. We shall find in the second 
tier of this study that the large event counts required 
to drive a fitting of the LSP mass will instead naturally 
favor the ULTRA selections. 

The missing transverse energy which has been dis- 
cussed is defined formally as follows. 
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Our simulations indicate a strongly peaked Poisson-like 
distribution of the histogram on missing transverse en- 
ergy, similar to the power spectrum of black body ra- 
diation with various temperatures. Figure Set ^ un- 
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FIG. 1: We depict the statistical distribution of missing energy events, per 20 GeV bin, scaled down to 1 fb~ for LHC center- 
of-mass collision energies of ^/s = (7, 8, 10, 12 and 14) TeV. Each frame displays the histogram count for 17 representative 
J--SU{5) benchmark points. The properties of the each individual point are listed in Table HI Also included for comparison 
is the "Snowmass Points and Slopes" (SPS) benchmark point SPSS 54]. The thick blue histogram peaking at low energies 
represents the ti + jets background sample. The lower right frame identifies the peak event count as a function of //™'°° for 
a continuous progression of LHC center-of-mass collision energies. The legend associates the gaugino mass M1/2, vector mass 
Mv, and LSP mass rriLSP of each of the 17 J--SU{5) points and SPS3 point with its respective color in the diagrams. 



mistakably highlights this physicahy analogous behavior 
by the way of five histograms delineating the > 9 jets 
event rate as a function of iJ™"^'', for LHC center-of-mass 
collision energies of ^/s = (7, 8, 10, 12 and 14) TeV, for 
each of the 17 representative J^-SU{5) points and the 
"Snowmass Points and Slopes" (SPS) benchmark point 
SPSS [5Jj . Curiously, the central value of this distribu- 
tion appears to correlate strongly with a small multiple 
of the LSP neutralino mass, a simple observation which 
shall develop into the cornerstone of our second tier anal- 
ysis strategy, to be presented subsequently. 

AU 2-body SUSY processes are included in Monte 
Carlo collider-detector simulation, which has been per- 
formed using the MadGraph [sl, [5a| suite, including the 
standard MadEvent ^, PYTHIA " and PGS4 ^ chain. 
Post-processing has been performed by a script CutLHCD 
of our own design (available for download (60|) that im- 
plements the desired cuts, and counts and compiles the 
associated net statistics. We likewise offer for down- 



load [6l[ the pair of Mathematica notebooks used to per- 
form the primary LSP analysis and plot generation for 
the present paper. The missing energy cut is relaxed for 
purposes of plotting histograms in _ff™"^'', and in particu- 
lar for the second tier analysis to follow. It is maintained, 
however, for the purposes of our first tier signal discover 
analysis. Our SUSY particle mass calculations have been 
performed using MicrDMEGAs 2. 1 [62|, employing a pro- 
prietary modification of the SuSpect 2.34 [63| codebase 
to run the RGEs. 



B. The SUSY CDM Discovery Index 

A well known indicator of the statistical significance 
of a collider signal in comparison to some competing 
background is the ratio S/\/B + 1 of signal events S to 
the square root of background events B, plus one. The 
"plus one" limiter in the denominator is a sanity check 
in the case of extremely small backgrounds, which pre- 
vents the signal count S from being numerically smaller 
than the target ratio. The conventional wisdom is that 
S/\/B + 1 > 5 should generally be considered favorable, 
with likewise a minimum of five signal events. Scaling 
both S and B by the factor N, the quadratic equation 
may be readily inverted to establish the following rela- 
tionship. 
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In the large-background limit, where the "plus one" term 
is not required, this reduces to simply 25B/S^. In Fig- 
ure ([2]), we have graphically extrapolated the necessary 
scaling N of integrated luminosity which would be nec- 
essary to achieve the target ratio of 5 for our discovery 
index. Since we have a well defined background count 
only for the current -^i = 7 TeV beam, we have ex- 
trapolated the observed event scaling for our own T- 
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FIG. 2: The projected minimal beam luminosity required 
to achieve a ratio S/\/B + 1 = 5 of signal to background 
events with at least nine jets and at least 150 GeV of miss- 
ing transverse energy is plotted for collision energies of y^s — 
(7, 8, 10, 12 and 14) TeV, for each of the 17 representative 
J--SU{5) benchmark points. Also included for comparison, 
with the diamond data markers, is the "Snowmass Points and 
Slopes" (SPS) benchmark point SPSS HI]. 



SU{5) Monte Carlo onto the backgrounds. Averaging 
across the model space, and omitting the heaviest two 
benchmarks (which scale somewhat more sharply than 
the primary spectra, demanding larger energies for any 
substantial event production), we find event ratios of ap- 
proximately (1.0 : 2.3 : 8.3 : 21 : 41) for beam energies of 
y/s = (7, 8, 10, 12, 14) TeV. As demonstrated in the plot, 
the logarithm of the scaling factor varies roughly linearly 
with the LSP mass, translating even slight differences 
into dramatic distinctions in the prospects for discovery. 
We likewise omit the heaviest two J'-SU{5) benchmarks, 
as well as the mSUGRA control point, from the linearized 
fitting. Incidentally, some portion of the visible deviation 
from logarithmic linearity in the figure may be attributed 
to the "plus one" limiter. 

To summarize the findings depicted in Figure ([2]), 
1 fb^^ of data at -^s = 7 TeV could be sufficient to 
unearth an LSP mass up to about 80 GeV, while a 
y/s = 14 TeV beam approximately doubles the discov- 
ery threshold to around 160 GeV, encompassing almost 
the entirety of No-Scale J^-SU{5), for the same lumi- 
nosity. This is consistent with our conclusions, pub- 
lished elsewhere [Ig, [23], that the model space below 
Ml/2 - 480 GeV and ttilsp - 90 GeV is disfavored 
by the recent CMS [11] and ATLAS [li] reports on 
the first femtobarn-plus of data, while the models near 
Ml/2 — 520 and ttilsp — 100 GeV are not only allowed, 
but well described by the low-statistics excesses which 
have been observed. Likewise, it is expected that the al- 
ready collected 5 fb~^ of data should be sufficient for con- 



elusive validation or refutation of that same model region. 
For comparison, the discovery threshold for the SPSS 
benchmark is overlaid onto the same plot, with diamond 
shaped markers. The prospects in this case are much 
less favorable than for the leading T-SU{5) benchmarks, 
with a massive target luminosity of more than 5, 000 fb~ 
at y/s = 7 TeV suggested for discovery in this case. The 
primary reason for this is simply that the LSP mass of the 
SPSS example is quite heavy, at more than 160 GeV. This 
is not, however, incidental. In fact, it is fairly generic that 
experimentally viable mSUGRA/CMSSM parameteriza- 
tions require rather heavy LSP particles in conjunction 
with a mass spectrum which is sufhciently globally heavy 
to evade the limits on squark masses. This in turn makes 
the discovery of a missing energy collider signal substan- 
tially less likely. On top of this, the mSUGRA/CMSSM 
models generally underperform the No-Scale T-SU{5) 
models in the ultra-high jet multiplicity regime, although 
it is difficult to make an absolute comparison, as the 
detailed spectra may be hierarchically distinct, and not 
only different in overall scale. In particular, we note that 
the example spectrum of Table [Tl] should not be naively 
compared against squark limits which are based on an 
assumed CMSSM/mSUGRA construction; a preferable 
approach is the direct fitting of a Monte Carlo simula- 
tion against detailed collider measurements which feature 
identically applied selection cuts [16, 29]. 

It is unsurprising to note that detectability is improved 
dramatically as one looks beyond the current ^/s = 7 TeV 
LHC to the future y/s = 14 TeV. Of course, the alternate 
path to achieving enhanced detectability, namely light- 
ness of the LSP mass, is one whose fundamental mecha- 
nism is outside of our hands. It is, however, precisely the 
combination of a light LSP with heavy squarks (excepting 
the light stop), again the central feature of the charac- 
teristic No-Scale T-SU{5) mi_^ < rrig < rrig mass hier- 
archy, which lies at the heart of this particular model's 
best prospects for imminent discoverability. The same 
model characteristic is simultaneously responsible for the 
preponderance of ultra-high multiplicity jet events which 
may further facilitate that discovery. Given the relatively 
small target luminosities required by the intermediately 
heavy model space, there is moreover a tangible possi- 
bility that the first delicate bouquet of a SUSY signal 
already perfumes the published data [16[, or lacking ap- 
propriate selections, that it remains enveloped beneath 
the detritus of lower jet count SM events [29] . 

At a time when the LHC, LEP, Tevatron and 
XENONIOO results are strongly impinging upon the key 
parameter space of the CMSSM, and many models which 
are favorable in terms LSP production are being ruled 
out by their squark spectra (which are likewise relatively 
light across the board), the efficient simultaneous i) eva- 
sion of experimental limits and ii) preservation of robust 
near term term testability represents no mean feat. Al- 
though we have focused here predominantly on collider- 
detector methodology, one should not meanwhile neglect 
the very real possibility that the direct detection exper- 



iments such as XENONIOO [l7| may trump the LHC to 
discovery and/or classification of a SUSY CDM candi- 
date [l3,[23- We have remarked [23], in particular, that 
the dual designation of Supersymmetric Dark Matter in 
fact requires the contribution of both search approaches, 
simultaneously complements and competitors. 



III. TIER 2: SUSY CDM CLASSIFICATION 
A. The Boltzmann Analogy 

The second tier of our analysis looks beyond discov- 
ery to model classification, and is focused on the detailed 
spectroscopic distribution of the missing energy counts. 
We were struck upon initial observation of the histogram 
counts against missing transverse energy i/^"**^, as dis- 
played in Figure Set ([T]), by the strong resemblance to 
the curves describing the power spectrum of a black body 
radiator at various temperatures. The memory stirs by 
this analogy to recall the treasure of information enci- 
phered in the cosmic microwave background which was 
extracted by application of this single same key. We have 
been most interested to further note the rough visible cor- 
relation between the distribution peak (i/™'^^) , and 
a small multiple of the LSP mass. Given the allowed 
kinetic excess over the jet invariant mass, the trigono- 
metric reduction from extraction of the beam-transverse 
component, the certainty of partial cancellation between 
the dual neutralino signal (with randomly oriented direc- 
tionality due to the multi level decay cascade), and the 
possibility of some substantial detector mismeasurement 
in at least one of the > 9 event jets, it does not at all seem 
obvious that any such clear correlation should persist. 

If one hopes nonetheless to reverse the correlation 
which is indeed observed, extrapolating from the statis- 
tics accrued under actual collider-detector operational 
conditions to the active LSP mass, then it is essential 
to have a quantitative understanding of the parameter 
dependencies in the event count function. An empirical 
reconstruction from the Monte-Carlo may serve as well as 
an analytically established relationship, if one is careful 
to choose a physically motivated ansatz. We will indeed 
here establish just such an empirical fitting relationship, 
led by the Boltzmann analogy to consider a distribution 
of the Poisson type. This single function will be demon- 
strated to simultaneously describe the entirety of the T- 
SU{5) model space. While the fitting function is not 
truly universal, insomuch as it fails to adequately model 
the SPSS control sample, we interpret this as an opportu- 
nity to distinguish experimentally between the signatures 
of various candidate scenarios. While the specific numer- 
ics may be quite model dependent, we speculate that the 
underlying Poisson correlation between the dark matter 
particle mass and missing transverse energy may be fully 
general. 

Most notably, we will successfully establish a specific 
numerical linkage between the LSP mass and the value 



of (if™'**^) , with cumulative uncertainties totaling ap- 
proximately 12 — 15%, which we find to be quite satisfac- 
tory. The co-linear relationship connecting the universal 
gaugino mass M112 with the LSP neutralino mass in No- 
Scale F-SU{b) must also be stressed, to the extent that 
precision measurement of the LSP mass via a global fit 
of the H!f^^^ histogram is virtually equivalent to an ex- 
perimental determination of M1/2 itself. As mentioned 
in Section III Al the larger overall event multiplicities as- 
sociated with the ULTRA selection cut procedure provide 
a much greater data density and uniformity than can be 
achieved with the CMS HYBRID cuts, and are thus prefer- 
able for use in our second tier study. Again, the selection 
cut on H^^^^ is relaxed for the purposes of generating a 
distribution in this parameter. 



B. An Empirical Poisson Fit 

A black body, in the Boltzmann limit, is described by 
a rescaling of the Poisson distribution 



leading coefficient A carries the unit [1 fb/20 GeV]. 
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which represents an expansion of the binomial distribu- 
tion in the limit of rare events, giving the likelihood of 
making N observations against an expected baseline of /i 
observations. It is important in this context, however, to 
note that the independent variable is taken as fj, rather 
than N, and the constant factorial denominator may then 
be replaced by a more general leading constant A (not 
to be confused with the trilinear soft terms). The differ- 
ential black body intensity per frequency step is then 
described by (/i — >■ i'/T,N — >■ 3), or alternatively by 
{fi — >■ l/AT, A^ — 7^ 5) with respect to differential steps 
in the wavelength. 

By trial, an excellent fit may be obtained for 
the frequency-style distribution, while the inverse 
wavelength-style distribution seems much less well suited 
to model the simulation. However, the non-zero event 
count in the limit H!^^^^^ — > forces one to consider (at 
least) a linear transformation of the horizontal axis, to 
fj, = m X {H!f^^^/100 GeV) -I- b, with dimensionless slope 
and intercept parameters. Such a relinearization of the 
Poisson distribution is familiar to high-energy experimen- 
talists, for example, as applied to the detection probabil- 
ity for a horizontal axis representing some (potentially 
unknown) event cross section, with integrated beam lu- 
minosity as the slope parameter, and the background pro- 
cess count representing the intercept. 

The exponent A'^ = 3 is inessential to our purpose, and 
a sequential consideration of somewhat larger integral 
values provided in several cases a comparable fitting qual- 
ity; nevertheless, with no single strongly favored alterna- 
tive, we will opt to respect the physical analogy by reten- 
tion of the cubic power. Together then, we adopt a three 
parameter fitting function, written as follows, where the 
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Explicitly differentiating the prior function, and setting 
the result to zero, one readily solves for the location of 
the function maximum on the i/™'^** axis. Secondarily 
reinserting (^miss-jmax i^g^^j^ jjj^q gq^ ^^ a,n expression 

for the function maximum is likewise obtained. 
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Integrating across H^^^^^ , one may establish the net event 
count per inverse femtobarn of luminosity, a measure 
which corresponds to the overall model discoverability, 
as emphasized in the first tier of the present work. The 
result is expressed via the partial Euler Gamma function, 
or may be polynomial expanded, as follows. 
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A root-mean-square quantification CT^miaa of the Poisson 
function width may also be integrated in closed form. 
The result will simplify considerably once a numerical 
value is established for the intercept parameter b. 
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Best fits, using a least-squares deviation metric, 
were composed for the parameters (A, rn, 6), for each 
of the complete Monte-Carlo event simulations of our 
seventeen J^-SU{5) benchmarks, repeated for each 
of the LHC center-of-mass collision energies y/s = 
(7, 8, 10, 12, 14) TeV. Figure i^ depicts a representative 
sample of six benchmarks, each at ^/s — 7 TeV, demon- 
strating the broad suitability of this function. We do 
find that beyond approximately Afi/2 = 650 GeV, the 
low event cross-sections result in large detection fluctu- 
ations, rendering a fit against the quantity of collected 
statistics somewhat tenuous. We have therefore omitted 
the two benchmarks in this category from the present 
analysis. 

Given data for the three parameters of Eq. ([7]), across 
our full model space and for various collider energies, 
the next task in the development of a suitably broad 
empirical function is the secondary functional correlation 
of each of these parameters with the model and collider 



10 



Mu-. = 3S5 GeV, M^- = 3575 GeV 



M,p_ = 410 GeV, My = 925 GeV 





100 120 

/wlsp in GeV 



FIG. 4: We depict the secondary linear fit of the slope pa- 
rameter m of each benchmark J-"-Sf/(5) model as a function 
of the LSP mass for each considered collider beam energy. 



FIG. 3: The counts of missing energy events, per 20 GeV 
bin, scaled down to 1 inverse femtobarn, are shown as blue 
dots from our Monte Carlo simulation. Six representative 
benchmarks are selected for display, each at a collider beam 
energy of ^/s = 7 TeV. The solid blue line is the best three 
parameter Poisson fit to each isolated model. 



inputs. We turn first to the main slope parameter m. 
Figure ^ demonstrates the strong linear correlation of 
m with the LSP mass for fixed beam energy ^/s. In fact, 
the basic beam energy independence of the parameter m 
is exhibited by the concurrency of the five overlaid color 
coded sample sets. 

An equally satisfactory fit for m is possible against 
Ml/2, which is unsurprising as one expects the basic pro- 
portionality rriLSP c< Ml 12, but since our purpose here is 
evocation of the LSP mass, we prefer the initially stated 
parameterization. A weak secondary, though also appar- 
ently systematic, dependence of m is observed against the 
vector-like mass My and the top quark mass mt- How- 
ever, the statistical scatter is large, and we judge the 
additional complication to be unjustified by the meager 
improvement in overall fit quality. Given each of the prior 
facts, we settle on a final numerical fitting for the main 
slope parameter m which is given by averaging across the 
beam energy -^/s for the five least-squares linearizations 
depicted in Figure (j4]). 
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The quality of this fit may be characterized in a stan- 
dard manner by the proportional deviation measure ct, 
which we notate with an overline to distinguish the rela- 



tive scaling procedure, and define as shown following. 
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We find a proportional deviation in this case of 'Wm — 
0.043, or crudely speaking, about 4% disagreement be- 
tween the single model fits for the slope parameter m 
and the global numerical fitting for all models. 

Moving on to the main intercept parameter h, Fig- 
ure ([5]) demonstrates substantially reduced linearity 
against the LSP mass. The situation is not improved 
for alternate potential independent fitting variables. Al- 
though there is a recognizable systematic trend, there is 
a wholesale abandonment of the trend for certain sample 
points, and the overall scatter remains large relative to 
the rise and fall of the best fit slopes. We therefore opt in 
this case for a simple average, as represented by the con- 
stant horizontal lines in the graphic. Once again, there is 
no discernible beam energy dependence, and the average 
of averages yields the following dimensionless best fit. 



5 = 0.721 



(14) 



Despite the overall appearance of statistical scatter, the 
clustering around the average is still sufficiently tight as a 
relative measure to maintain a proportional deviation of 
only CTf, = 0.088, representing disagreement on the level of 
about 9% between the single model fits for the intercept 
parameter b and the global all model average. 

Finally, we take up the secondary fit of the main scale 
parameter A of Eq. ([7]), which is critical to the over- 
all model discoverability. We observe, as attested to 
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FIG. 5: We depict the secondary averaging of the intercept 
parameter b of each benchmark J^-SU{5) model as a constant 
function of the LSP mass for each considered collider beam 
energy. 



by the first plot of Figure Set ^ that it is the natu- 
ral logarithm of A x (1 GeV/1 fb) which scales linearly 
with the LSP mass, rather than A itself, implying an 
exponential dependence on the dimensionless parameter 
(m-LSp/100 GeV). A fit of comparable quality was ob- 
tained in a log-log fitting procedure, but the asymptotics 
of the power law scaling were deemed less physical. The 
second plot of Figure Set dH) , with the A axis linearized, 
graphically demonstrates the incredibly large dynamic 
range of the A parameter across our model and LHC 
beam energy space, which spans a full three orders of 
magnitude. It is apparent additionally in both plots that 
there is a tertiary functional dependence of the fit for A 
against the LSP mass on the beam energy y^, particu- 
larly in the fit intercept, but also in the slope. This is 
to be expected, of course, insomuch as detection counts 
should scale rapidly with the increasingly large interac- 
tion cross sections afforded by the increased phase space 
of a more energetic collider, as depicted in the lower right 
frame of Figure Set ([1]), which identifies the peak event 
count as a function of H!f^^^ for a continuous progression 
of LHC center-of-mass collision energies. 

The plots of Figure Set ([6]) demonstrate moreover a 
substantial spreading, within a single beam energy con- 
tour, for various benchmarks of similar LSP mass. This 
fact, together with the extreme exponential scaling, sug- 
gest that if a suitable fitting is possible, it must also take 
into consideration other input parameters such as Aiy, 
tan/3, and mt. We omit M1/2 from this list because of 
its strong co- linearity with rriLSP- Since the distribution 
peak and width are unaffected by the overall scale term 
A, we are simultaneously less reluctant in this case to ad- 



mit a somewhat more intricate parameterization. Guided 
by several additional plotting test cases for subordinate 
dependencies (the printing of which we here suppress), 
a persistent eye for maximal functional simplicity and 
physicality, plus no small measure of trial and error, we 
settle at last on the following form for the numerical value 
of A, in [1 fb/20 GeV] units. 
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The global fit against the main scale parameter A is mea- 
sured to have an extraordinarily small proportional de- 
viation of a A = 0.036, i.e. less than 4%. With this fi- 
nal piece in place, the quality of the overall fitting for- 
mula in (tolsPi y/s, My, tan /3 and mt) which is expressed 
by the union of Eqs. ([71 [TH [14] and [15]) may be judged 
by visual inspection of Figure Set ([7]). Beam energies 
of ^/s = (7, 8, 10, 12 and 14) TeV are there superim- 
posed for each of four representative benchmark sam- 
ples. The blue dots represent the histogram binning from 
the Monte Carlo, the blue curves represent the least- 
squares fit to the single sample space in isolation, and the 
red curves represent the global fitting function. Figure 
Set ([8]) next depicts a scan across all benchmarks (except- 
ing A/1/2 > 650 GeV) for center of mass beam energies 
of 7i = 7 TeV and ^/s = U TeV. The histogram data 
and individual sample fits are omitted for visual clarity. 
The blue rings denote the distribution peak. 



C. The Once and Future LHC 

Having established the global Poisson fitting function, 
we now pause to reflect on the key purpose of this sec- 
ond tier study, granting for argument's sake an immi- 
nent SUSY collider signal observation, and looking by 
the mind's eye a bit farther ahead, perhaps beyond the 
present LHC to a forthcoming -y/s = 14 TeV incarnation, 
and a time after the initial excitement of the SUSY CDM 
discovery has given way to the business of extracting the 
distinguishing details of the physical SUSY projection. 

We consider first the parameter subset of just m and 
6, as they are alone sufficient to fully characterize the 
location of the Poisson peak in iJ™"''' from Eq. ([8]), as 
well as the distribution width of Eq. (jlip . Substituting 
back into each of these relations, we obtain the follow- 
ing, which are expected to apply across the full T-SU{5) 
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FIG. 6: We depict the secondary exponential fit of the slope parameter m of each benchmark T-SU{5) model as a function of 
the LSP mass for each considered collider beam energy. The left-hand plot is logarithmic, while the right-hand plot is linear. 



model space. 
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(2.04 







(16) 
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The expression from Eq. (ITU)) for the net event count 
per inverse femtobarn of luminosity similarly reduces as 
shown following, where (m. A) carry the functional de- 
pendencies outlined in Eqs. (|12I15I) . 



/ Events \ „„ „ /A 
r ^ 29.8 X — 



(18) 



The overall proportional deviation from the actual Monte 
Carlo counts is less than 1%, using the local {A,m,b) 
parameters, and still less than 5% using the global fitting 
parameters. 

Inverting Eq. (|16p. we can get the desired expression 
for the LSP mass in terms of the observed location of 
the Poisson peak, or alternatively along with Eq. ([S]), in 
terms of the Poisson fitting parameters. 



TOLSP ^ 204 GeV - 163 GeV x 



= 204 GeV - 163 GeV x 



/ 100 GeV 



3-6 



(19) 



Comparing our input LSP masses, as processed through 
the collider-detector Monte Carlo simulation, to the value 
obtained by fitting the resulting H^^^^ distribution ac- 
cording to the Eq. ([T]) ansatz, and passing the fit param- 
eters through Eq. (|19p . we observe a proportional devi- 
ation of CTrriLSP = 0.058, or about 6% error, induced by 



the closed circuit. In a realistic application, this would of 
course necessarily be compounded with any uncertainty 
in the fitting procedure itself. 

Uncertainty in the overall fit against data comes from 
multiple sources. First, there are statistical fluctuations 
and detector errors in the measurements of the i/™'^** 
histogram itself. These should be reasonably modeled 
in our approach, which again includes a post processing 
phase by PGS4 [53 to simulate the necessarily faulty and 
incomplete nature of real world data. However, it should 
be noted that although our histograms are normalized to 
1 fb~ of data, we have generally overrun the simulation 
by some (possibly large) factor and scaled down, which 
will tend to reduce random scatter. Secondly, it is not a 
priori obvious that Eq. ([7]) itself carries any fundamen- 
tal validity, although the fitting exemplified in Figure ^ 
does speak strongly for its use. In an attempt to quan- 
tify these two potential sources of error, we computed, for 
each T-SU{5) model and beam energy combination, the 
absolute standard deviation between the Monte Carlo bin 
heights and the best fit Poisson curve of the given model, 
dividing at the end by the average bin height, and fi- 
nally averaging across the full model space. An error of 
about 12% was observed. Using instead the global fit of 
Eqs. (|12ll4p to generate the (m, b) parameters based on 
the known LSP mass, the error increased to about 15%. 

An additional, and potentially more problematic, 
source of real world error, which we have easily 
sidestepped in simulations, is the difficulty of deconvolv- 
ing the SM background distribution from the SUSY sig- 
nal. We argue in Refs. [2J, ^^ that the substantially 
elevated jet cut of greater than or equal to nine jets over- 
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niLSP = 75 GeV, «„, = 395 GeV. M,. = 2075 GeV. tiiii(/!) = 19.7, m, = 172.5 GeV 



«:ls|i = 82.8 GeV, «„, = 445 GeV, «,. = 1 125 GeV, tiin(/J) = 19 



', = 174.4 GeV 





mLSP = 100,5 GeV, M,,, = 505 GeV, M,- = 3700 GeV, Hn(^) = 21, m, = 172,6 GeV 



mLSP = 133,2 GeV, M,,, = 650 GeV, My = 4700 GeV, Hn(^) = 22, m, = 173,4 GeV 





FIG. 7; The counts of missing energy events, per 20 GeV bin, scaled down to 1 inverse femtobarn, are shown as blue dots 
from our Monte Carlo simulation. Four representative J--SU{5) benchmarks are selected for display, for collider beam energies 
of v^ — (7, 8, 10, 12 and 14) TeV, from bottom to top. The solid blue lines are the best three parameter Poisson fits to each 
isolated model. The solid red lines are the global fit based only on the input parameters (niLSP, \/s, Mv,tan/3 and nit). 



whelmingly compromises the observability of the SM sig- 
nal. However, the footprint of the new physics is itself 
sufhciently narrow that our simulation of the the it + jets 
background (in addition to the QCD and vector boson 
plus jets components which we have not personally mod- 
eled, but which are accounted by the large collaborations 
such as CMS) can generate a sufficient number of events 
to encroach upon the J^-SU{5) signal. As exhibited by 
the ti distribution in Figure Set ([J), the SM background 
contributions may be expected to peak at lower energies, 
and a hard lower bound cut on H^^^^ is thus quite effec- 
tive if one is interested in the integrated signal, for the 
pure initial purposes of a proof of signal existence. In 
the present case, however, it is quite troublesome that 
the missing energy (//miss-jmax pgg^]^ should occur in the 
vicinity of 100 GeV, potentially shrouded and obscured 
by the SM veil. 

However, there is some substantial comfort available 
from the explicit analytic functional form of Eq. (I7|) . If 
one tacitly accepts this ansatz distribution, then the be- 
haviors at the peak and distance from the peak are im- 
plicitly coupled. In fact, it is possible to fully extrapolate 
the precise form of the three parameter fit in the region 



where background competes strongly only by carefully 
observing the distribution in the region with low back- 
ground. Such an approach, combined certainly with a 
careful subtraction of the estimated SM background at 
smaller iJ™"''' values, may present a surprisingly direct 
and well resolved window on the mass of the underlying 
LSP particle which is responsible for the missing energy 
signal. The essential decoupling of the simple Eq. P^ 
result from the collider beam energy seems to us quite 
remarkable. It is also worth remarking that the Eq. ([15]) 
expression of the scale parameter A, while not so for- 
mally simple as the expression for the slope parameter 
m, does contain corroborating information on the LSP 
mass, which is in some sense orthogonal to that extracted 
in Eq. (IE 



D. The General and The Proprietary 

In Figure Set ([9]), we present the parallel analysis of 
our SPSS mSUGRA control sample. The actual his- 
togrammed events (blue dots), along with the blue line 
best fit, are contrasted with the corresponding global J^- 
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niLsp = 161-7 GeV, M„, = 400 GeV, M„ = 90 GeV, lan(;3) = 10, ni, = 175 GeV 
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FIG. 8: The global fit to missing energy is shown as a solid 
red line for each J^-SU{5) benchmark model, with heavier 
spectra at the bottom, and lighter toward the top. The blue 
rings denote the distribution peak. The plots are based only 
on the input parameters (771lsp, \/s, Afv, tan/3 and rrit), and 
not directly on the Monte Carlo. Collider beam energies of 
y/s = 7 TeV and y/s = 14 TeV are shown for the sake of 
comparison. 



SU{5) fit, using the SPSS parameters, vifith A/y = 0. 
The divergence is clear, in both scale and location of the 
distribution peak, immediately demonstrating that the 
empirical fit does not directly generalize into alternative 
model spaces. Nevertheless, the essential Poisson charac- 
teristic appears to be persistent, and we speculate that 
this bulk feature may be broadly applicable, while the 
specific numerical coefficients may be model dependent. 
There is also a second sense in which the "tier 2" pro- 
cedure may be general in principle, though again propri- 
etary in practice. This is with respect to the prospects for 
overcoming competing backgrounds. The ultra-high jet 
multiplicity cutting procedure, which we suggest is capa- 
ble of sufficiently elevating the No-Scale JF-SU{5) H^''''' 
distribution above the SM background that a function 
fit might be possible, may not be as effective for generic 
mSUGRA/CMSSM samples, which may moreover suffer 
from a necessarily heavier LSP candidate. We must also 
here reemphasize two key points: (1) Measuring the LSP 
neutralino mass is essentially equivalent to measuring the 




Vj"= 10 TeV 
mLsp = 161.7 GeV, M„,= 400 GeV. M„ = 90 GeV, lanOT= 10, m, = 175 GeV 




FIG. 9: Statistical distributions of missing energy events, per 
20 GeV bin, scaled down to 1 fb~^ for LHC center-of-mass 
collision energies of ^/s = (7 and 14) TeV, for the mSUGRA 
benchmark SPSS. Each frame displays the histogram count 
as blue dots from our Monte Carlo simulation. The solid blue 
line is the best three parameter Poisson fit for this isolated 
model. The global T-SU{5) fit to missing energy is shown as 
a solid red line for the corresponding parameters of the SPS3 
benchmark, with Afv = 0. 



universal gaugino mass M1/2 in No-Scale J^-SU{5) be- 
cause the overall mass rescaling of the model is based 
only on this single dimensionful GUT scale parameter; 
(2) If our method for the LSP neutralino mass deter- 
mination cannot be broadly applied outside of No-Scale 
J^-SU{5), the absence of this rescahng property may be 
the deep reason. 

One should be cautious in attempting to extend the 
discrepancy manifest by the global fitting function in 
Figure Set ©. For example, the overage of the red 
line fit, which might casually be interpreted to imply a 
generically weaker discovery potential for the mSUGRA 
candidates, can easily be converted to a severe under- 
estimate by the transformation of tan/3 from 10 to 30, 
a switch which roughly corresponds to adoption of the 
SPSlb Snowmass benchmark, rather than SPSS. In this 
example, the strong dependence of Eq. (fT5|) on the ratio 
tan/3, which is softened in J^-SU{5) by a generic prox- 
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imity to the numerical value of 20, is exaggerated for 
the SPS models. The lesson here is that comparisons 
of this sort are subtle, and that the conventional wis- 
dom developed in one model space does not necessarily 
transfer intact to an alternate class of constructions. We 
remark, however, that the clear model dependency of the 
applicable Poisson fitting function, itself contains essen- 
tial information, which may ultimately prove helpful in 
the process of distinguishing between contending efforts 
to explain the ongoing accumulation of LHC data. 

As a further control on the cross-applicability of our re- 
sults, we fully regenerated the primary LSP fitting anal- 
ysis presented in Section IIIIBI for two related but dis- 
tinct scenarios, i) substituting the CMS HYBRID cutting 
methodology, and ii) retaining the ULTRA selections, but 
substituting a kt jet clustering algorithm (specifying an 
angular distance parameter of 0.4) for our default cone 
algorithm (these being the two analysis options avail- 
able within the PGS4 [531 detector simulation) . While the 
event count scale, as characterized by our parameter A, 
is affected in both cases, we find the correlation between 
tolsp and the maximum of the the event histogram in 
H^^^^, as characterized by our parameters m and 6, to 
be rather stable. Moreover, the leading modification to 
the scale parameter A is a constant multiplicative factor. 
This is a rather reassuring result, which preserves the 
correlation established by Eq. ([T5|. which is again like- 
wise independent of the collider beam energy. We find 
it quite satisfactory that the overall event scale, which is 
the most difficult parameter to fit, and the most sensitive 
to incidental variation of the experimental technique, is 
simultaneously the least critical carrier of actual informa- 
tion on the J^-SU{5) LSP particle and SUSY spectrum. 



IV. CONCLUSIONS 

The LHC era is rapidly evolving, with each major ex- 
perimental collaboration having already accumulated a 
reported 5 fb~ of integrated luminosity at the close of 
2011. Looking into the not so distant future, we must 
plan to focus not only on initial attainment of the much 
anticipated SUSY signal discovery, but also on mining 
the wealth of information present in the data for clues 
to the reconstruction of a detailed SUSY mass spectrum. 
We have concentrated in this presentation on precisely 
such a two tiered objective. Firstly, we have comprehen- 
sively documented the prospects for discovery of the full 
No-Scale J^-SU{5) model space at various LHC center- 
of-mass energies, based on the raw count of events with 
greater than 100 GeV of missing transverse energy and 
at least nine jets. Secondly, upon observing the dual cu- 
riosities of i) a material correlation between the i/™'^** 
histogram peak and a small multiple of the LSP mass, 
and ii) the strikingly familiar shape of the missing en- 
ergy distribution, reminiscent of the power spectrum of 
a black body radiator at various temperatures, we have 
implemented a broad empirical fit of our extensive Monte 



Carlo H^^^^ simulation against a Poisson distribution 
ansatz. We advanced the resulting fit as a theoretical 
blueprint for deducing the mass of the LSP, utilizing only 
the missing transverse energy in a statistical distribution 
of > 9 jet events. We judge the error directly associated 
with this technique to be within the satisfactory range 
of 12 — 15%. It will be interesting to see whether this 
supersymmetry mass extraction instrument might find 
some value by inclusion in the experimentalist toolbag, 
as we journey into and beyond the presumed discovery 
phases of the LHC operation in the next two years. 

We have presented this analysis in the context of the 
No-Scale J^-SU{5) model, distinguished by its essential 
phenomenological consistency, profound predictive ca- 
pacity, singularly distinctive experimental signature, and 
imminent testability. The core of the model space around 
and above M1/2 = 500 GeV and mLSP = 100 GeV re- 
mains experimentally unblemished, in stark contrast to 
the incisive cuts already made into the heart of the more 
traditional mSUGRA/CMSSM formulations. This deft 
evasion of all existing LEP, Tevatron, and LHC sparticle 
mass limits stems from the characteristic sparticle mass 
hierarchy mj:_^ < rrig < mq of a light stop and gluino, 
both comfortably lighter than all other squarks, a feature 
which is stable across the full model space. With some 
reports claiming that less than one percent of the for- 
merly viable mSUGRA and CMSSM parameter space has 
withstood even the earliest accumulation of LHC data, 
the existence of a clearly defined, deeply theoretically 
motivated model which simultaneously exhibits consis- 
tency with all current experiments alongside testability 
at near term future experiments, is of some pressing in- 
terest. The seventeen example spectra which we have 
studied span the model space of bare minimal constraints 
on No-Scale J^-SU{5), and establish a material correla- 
tion between the mass of the dark matter particle, taken 
to be the LSP neutralino, and the collider missing energy 
signal. It appears that the stability of the SUSY mass 
hierarchy, with the overall scale driven by a single dimen- 
sionful parameter at the GUT scale, may be responsible 
for facilitating the remarkably robust and simple fitting 
relationship which has been observed. Together, these 
facts speak efficiently to the suitability of No-Scale J^- 
SU{5) as a leading contender in the race to describe post 
Standard Model physics at the dawn of the LHC era. 
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